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COMPUTER APPARATUSES AND PROCESSES FOR ANALYZING A SYSTEM 

HAVING FALSE START EVENTS 

This invention was made under DOE CRADA No. LA95C01191. The 
Government has certain rights in the invention. 

CROSS REFERENCE TO RELATED APPLICATION 

This application claims benefit and priority to provisional application 60/202,010 
filed May 4, 2000. 

TECHNICAL FIELD OF THE INVENTION 

The present invention relates to the field of processes and apparatuses for 
analyzing a system, and, more particularly, to the field of processes and apparatuses for 
simulating and analyzing repairable systems, such as manufacturing systems, queuing 
systems, etc. 

BACKGROUND OF THE INVENTION 

Reliability is a measure of the probability that parts, components, products, or 
systems will perform their designed-for functions without failure in specified 
environments for desired periods at a given confidence level. Typically, reliability is 
expressed as a decimal fraction (e.g., 0.832). Reliability Engineering encompasses the 
theoretical and practical tools by which the probability and capability of parts, 
components, equipment, products, and systems to perform their required functions can be 
specified, predicted, tested, demonstrated, installed, and initialized. Powerful reliability 
engineering tools which enable the accurate prediction of the reliability of parts, 
components, and systems can provide a company with a significant competitive 
advantage. For example, accurate prediction of the reliability of a proposed or existing 
manufacturing or production line can decrease costs, increase speed to market of new 
products, and provide more predictable project outcomes. Accurate reliability prediction 
can also be used to identify and allocate resources for the implementation of process 
changes which can increase the reliability of a manufacturing system. 
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changes sequentially through time) to simulate the reliability of production systems or 
their subsystems. For use in simulating the dynamics of a complex manufacturing 
system, these simulations can be time consuming, and expensive to develop and run. 
Simplifications can introduce errors which challenge the validity and usefulness of these 
5 simulations. Benefits in speed to market, right to market, and predictable project 
outcomes, as above, may be had from methods which offer high degrees of accuracy in 
simulating the performance of existing systems and may be used for simulating scenarios 
which intend to alter such systems, or to simulate systems which may use new 
combinations of subsystems from valid models of existing manufacturing systems. 

10 Manufacturing costs and capacity can be better understood and controlled with more 
accurate prediction methods resulting in less market upset, particularly during the early 
phases of a new product introduction. 

As such, there is a need for simulation methods which can predict the coupling or 
compounding effects that premature failures, or "stop-and-go" operations have on the 

15 manufacturing system performance. Further, there is a need for simulation methods 
which enable better utilization of downtime for repair and restoration of system operation. 
That is, repair actions are taken specifically to raise the probability of a successful restart 
of the failed subsystem, and to reduce the frequency of false starts, and short times to 
failure due to inadequate repair (i.e., maintenance induced failure). 

20 

SUMMARY OF THE INVENTION 

Computer apparatuses and processes for analyzing a system are provided. The 
apparatuses and processes incorporate steps of collecting data from a first system, wherein 
25 the first system has a plurality of failure modes and the data relates to the failure modes, 
parameterizing the data for use in a simulation of a second system, and executing a 
simulation of the second system. The step of executing the simulation also includes 
determining whether the second system will encounter a first false start event based upon 
the data collected from the first system. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

While the specification concludes with claims particularly pointing out and 
distinctly claiming the invention, it is believed that the present invention will be better 
5 understood from the following description taken in conjunction with the accompanying 
drawings in which: 

Fig. 1 is a schematic illustration of exemplary process in accordance with one 
aspect of the present invention; 

Fig. 2 is a schematic illustration of an exemplary paper towel manufacturing 

10 system; 

Fig. 3 is a schematic illustration of an illustrative hierarchy for the manufacturing 
system of Fig. 2; 

Fig. 4 is a plot of exemplary downtime system loss event data; 

Fig. 5 is a schematic illustration of networked computers suitable for use with the 
1 5 present invention; 

Fig. 6 is a schematic illustration of a preferred architecture for a simulation made 
in accordance with the present invention; 

Figs. 7 to 9 are schematic illustrations of a preferred process implemented using 
the simulation of Fig. 6; and 
20 Fig. 10 is a schematic illustration of a preferred process implemented for the 

accumulator of the simulation of Fig. 6. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

25 Reference will now be made in detail to the present preferred embodiments of the 

invention, examples of which are illustrated in the accompanying drawings wherein like 
numerals indicate the same elements throughout the views. For purposes of clarity and 
simplicity, the following terms are used herein. 

As used herein, the term "system" is intended to refer to any set of components, 

30 processes, operations, or functions which deliver a product or service. Exemplary 
systems include manufacturing and production systems, airline flight operations, 
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maintenance operations, queuing operations (e.g., traffic control, reduction of waiting 
time in lines or queues such as bank teller windows, gas pump, grocery checkout, etc.), 
industrial and military asset deployment (e.g., troop deployments, air wing sizing, rental 
car and truck fleet sizing, warehouse and transportation logistics). 
5 As used herein, the term "loss event" is intended to refer to any event which 

negatively affects the performance or function of a system or one of its components (e.g., 
system stops or component downtime, product quality decreases, increased operator 
intervention, idle trucks or cars, empty lines or queues), wherein each loss event has a 
cause and a failure mode associated therewith. 

10 As used herein, the term "failure mode" is intended to refer to a description of 

how a system can fail to perform its intended function. Each failure mode can have 
associated therewith one or more causes; a downtime during which the failure mode has 
caused a system, or component thereof, to cease to perform its intended function; and an 
uptime until the next occurrence of a failure mode, during which time a system, or 

1 5 component thereof, is performing its intended function. 

As used herein, the term "cause" is intended to refer to the reason why a failure 
mode occurs. 

As used herein, the term "life" is intended to refer to the time from a starting point 
to the occurrence of an end point. 

20 As used herein, the term "competing cause failure mode" (CCFM) is intended to 

refer to a failure mode which occurs according to its own distribution, which is 
independent of the distributions of other failure modes and where the competing cause 
failure mode's uptime terminates at the occurrence of any loss event. In simulations of 
the present invention, the uptime of each competing cause failure mode is regenerated 

25 after the occurrence of any loss event as if anew. 

As used herein, the term "cumulative cause failure mode" (CMFM) is intended to 
refer to a failure mode which occurs according to its own distribution, which is 
independent of the distributions of other failure modes and where the cumulative cause 
failure mode's uptime is not affected by the occurrence of loss events unrelated to the 

30 subject cumulative cause failure mode. In simulations of the present invention, the 
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uptime of each cumulative cause failure mode is only regenerated anew after its 
expiration. 

As used herein, the term "false start event" is intended to refer to a loss event 
which occurs quickly relative to the expected life of a system (e.g., two minutes or less 
5 for a manufacturing system which may have an average expected life of twenty to thirty 
minutes) after a system, or component thereof, resumes acceptable operation. 

As used herein, the term "series system" is intended to refer to a system which 
fails if any of its components fail. 

As used herein, the term "parallel system" is intended to refer to a system which 
1 0 fails if all of its parallel components fail. 

As used herein, the term "uptime" is intended to refer to the period of time or life 
of a system, or component thereof, during which it is performing its intended function. 

As used herein, the term "downtime" is intended to refer to the period of time or 
life of a system, or component thereof, during which it is not performing its intended 
1 5 function due to a loss event. 

As used herein, the term "availability" is intended to refer to ratio of total uptime 
for a system to the sum of the total uptime for a system plus the total downtime for a 
system. 

As used herein, the term "Mean Time Between Failures" (MTBF) is intended to 
20 refer to as the ratio of the total uptime for a system over the total number of loss events 
for the system. 

As used herein, the term "Mean Time To Repair" (MTTR) is intended to refer to 
ratio of the total downtime for a system over the total number of loss events for the 
system. 

25 As used herein, the term "parameterize" is intended to refer to the process of 

characterizing or fitting data according to a parametric equation (e.g., equations 
containing parameters of shape, scale, and location to describe data, such as Weibull 
equations, log-normal equations, normal equations, etc.). 

Referring now to Fig. 1 and in accordance with one aspect of the present 

30 invention, a process for analyzing and/or simulating an exemplary system will now be 
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described. The process 20 begins at step 21 where a hierarchial level of the system is 
selected. The hierarchial level defines the system to be analyzed and represents the depth 
to which the components of a system are simulated. For ease of discussion, the present 
invention will be described herein with respect a system in the form of a manufacturing 
5 system and the machines which form part of the manufacturing system. More 
particularly, the present invention will be described herein with respect to a series towel 
manufacturing system 22 as schematically illustrated in Fig. 2. The manufacturing 
system 22 comprises a sheet roller 24 having wound thereabout a paper sheet 26 having a 
finite length. The sheet 26 passes through an embosser 28 which imprints a pattern upon 

10 the sheet 26. The sheet 26 is cut in the machine direction into discrete portions 32 at a 
first cutter 30. The portions 32 are wrapped about a cylindrical paper core at the core 
roller 34 after which the rolled portions 36 are deposited in an accumulator 38. The rolled 
portions 36 are next cut in the cross machine direction at the second cutter 40. The twice 
cut rolled portions 42 are delivered to a wrapper 44 which wraps a plurality of the rolled 

15 portions 42 in a polymeric film 46 which is drawn from a polymer roller 48. A plurality 
of sensors 50 are distributed throughout the manufacturing system 22 to detect loss events 
of the manufacturing system, as discussed more fully hereafter. While the manufacturing 
system 22 is described herein at a machine hierarchial level, it will be appreciated that the 
manufacturing system 22 can be described at higher and lower hierarchial levels based 

20 upon the analysis objective and desired level of simulation accuracy and that these 
various levels can be simulated collectively or individually as desired. Fig. 3 sets forth 
various illustrative hierarchies above and below that described for the manufacturing 
system 22 of Fig. 2. For instance, a higher hierarchial level might include a distribution 
network 52 comprising a plurality of manufacturing plants 54 (a parallel system), each of 

25 which have a plurality of the manufacturing systems 22, wherein the manufacturing 
systems 22 have a plurality of machines 56 (e.g., the roller 24) associated therewith and 
each machine 56 has a plurality of sub-components 58 (e.g., a roller bearing for the sheet 
roller 24). Further, while the present invention is described herein with respect to a 
manufacturing system, it will be readily appreciated that other products, services, 
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manufacturing systems, and systems are suitable for use with the present invention as 
previously described. 

Referring again to Fig. 1, the system hierarchy selected at step 21 is analyzed to 
identify its cumulative and competing failure modes and their causes. Set forth in Table 1 
5 below is an illustrative listing of cumulative and competing failure modes and their causes 
for the manufacturing system 22. 



TABLE 1 



Machine 


Failure 
Mode No. 


Failure Mode 


failure Mode 
Type 


Failure Mode Causes 


*io. of downtime 
distributions 


Sheet Roller 


1 


Sheet jam 


Competing 


Sheet break — hole in 
sheet 

Sheet break - glue 
buildup 


3 




2 


Roller change 


Cumulative 


Wear out - bearing 
failure 


1 




3 


Belt failure 


Cumulative 


Belt break 




Embosser 


4 


Infeed jam 


Competing 


Glue on roll 




First Cutter 


5 


Blade failure 


Competing 


Drive belt slipping 




Core Roller 


6 


Core 

misalignment 


Competing 


Short core 




Accumulator 


7 


Accumulator 
full 


Competing 


Accumulator full 




Second cutter 


8 


Blade failure 


Competing 


Blade wear 




Polymeric 
roller 


9 


Roller change 


Cumulative 


Roll expired 




Wrapper 


10 


Polymer sheet 
jam 


Competing 


Drive setup wrong 






11 


Polymer sheet 
alignment 


Competing 


Registration marks not 
aligned 





By way of example with respect to Table 1, the sheet roller 24 has three failure modes: 
10 two cumulative cause failure modes and one competing cause failure mode. Whether a 
failure mode is cumulative or competing can be based upon the cause of the failure mode 
or rules of construction (e.g., all failure modes are competing unless identified otherwise). 
The rules of construction can be tailored according to the objectives or use of the 
simulation. The sheet jam failure mode is competing, because, sheet instabilities can 
15 cause a sheet jam anew after repair of any loss event. The belt failure and roller change 
(i.e., replacement of the roll because the sheet is depleted) failure modes are cumulative, 
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because they are dependent upon the finite life of the belt and sheet irrespective of other 
unrelated loss events. The failure modes can have multiple causes (see, e.g., failure mode 
number 1 in Table 1) and downtime distributions. Further, the failure modes do not 
necessarily have a correspondence between the number of failure mode causes and the 
5 number of downtime distributions. The downtime distributions represent independent, 
identically distributed repair time characteristics for a single failure mode. For instance, 
the sheet jam failure mode has three downtime distributions, such as a short repair time 
period, a medium repair time period, and a long repair time period. Selection of the 
number of downtime distributions can be determined using methods known in the art, 

1 0 such as segmented regression wherein discontinuities between segments each represent a 
separate distribution. The probability C of falling within any given downtime distribution 
is determined by the ratio of the number of loss events within a downtime distribution 
over the total number of loss events for the subject failure mode. For example with 
reference to Fig. 4, an exemplary plot for the sheet jam competing cause failure mode of 

15 the sheet roller 24 is illustrated. The total number of loss events illustrated in Fig. 4 is 
188. For the first downtime distribution 70, the number of loss events is 109 and 
therefore the probability C of a downtime occurring in the first downtime distribution 70 
is 0.58. Similarly, the number of loss events for the second downtime distribution 73 is 
76 and therefore the probability C of a downtime occurring in the second downtime 

20 distribution is 0.40. 

Referring again to Fig. 1 and the process 20, time-based system loss event data 
(e.g., uptimes, downtimes) is collected from the manufacturing system 22 and 
parametrically analyzed at steps 60 and 62 once the hierarchial level of the manufacturing 
system 22 has been selected. The collection of this system loss event data is used to 

25 characterize and quantify the failure modes of the manufacturing system 22. The system 
loss event data which is collected preferably includes failure mode causes, time stamps 
marking the beginning and end of uptime and downtime periods, and the absolute time of 
uptime periods and downtime periods. Depending upon the failure mode cause, various 
detection techniques can be employed for collection of the data. Illustrative methods are 

30 set forth in Table 2 below. 
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TABLE 2 



Detection Method 


Description 


Sensor 50 equals failure mode cause 


The type and location of the sensor 50 is specific to a 
single failure mode cause 


Sensor 50 approximately equals the failure mode 
cause 


Multiple failure mode causes may trigger the sensor 
50, but generally one failure mode cause 
predominates 


Sensors 50 and logic/timing equals the failure mode 
cause 


Logic and timing combinations of sensors 50 can be 
used to deduce a single failure mode cause 


Operator input 


No reliable method based on sensors can detect the 
failure mode cause and an operator must make an 
entry 



As indicated, the system loss event data can be collected using the sensors 50 and/or 
programmable logic controllers (or PLCs). Illustrative system loss event data for the 
5 manufacturing system 22 is set forth below in Table 3. Enough data should be collected 
for each failure mode to establish a statistical sample. More preferably, at least about 3 
events are recorded for each failure mode. Most preferably, between about 1 5 events and 
about 30 events are recorded for each failure mode. The system loss event data can also 
include other information, such as a date stamp. While Table 3 is directed to system loss 
10 event data for a manufacturing system, other data types can be monitored for other 
systems, such as idle time for queue systems (e.g., gas pump available but idle). 

TABLE 3 



Event 


Uptime 
(minutes) 


Event beginning 
time stamp 
(minutes) 


Downtime 
(minutes) 


Event ending 
time stamp 
(minutes) 


Failure 
mode no. 


1 


25.1 


25.1 


6.8 


31.9 


6 


2 


35.4 


67.3 


5.1 


72.4 


8 


3 


36.8 


109.2 


2.5 


111.7 


1 


4 


14.3 


126.0 


3.2 


129.2 


11 


5 


0.4 


129.6 


10.3 


139.9 


6 


6 


1.4 


141.3 


8.6 


149.9 


2 



15 Once the system loss event data has been collected, it is next parameterized at step 

62 of the process 20. The system loss event data can be parameterized using one of 
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several statistical distribution/function types known in the art, such as an exponential 
distribution, a normal distribution, a Weibull distribution, or a log normal distribution. 
The Weibull distibution (in particular the Weibull Probability Density Function or PDF) 
is particularly useful for characterizing uptimes while downtimes are typically 
5 characterized according to a log normal distribution. False start events are preferably 
calculated according to the following equation: 

False start probability (P) = 1 -(total number of false start events for a failure mode) (1) 

total number of loss events for the system 

10 

Thus, a false start probability P of 0.97 (or 97%) means that 97% of the time, there is not 
a false start for this failure mode following any loss event for any failure mode. A typical 
Weibull PDF for uptime has the following form: 

15 Probability of Failure = (Beta/alpha) (uptime/Alpha) (Beta-1) e (-(((uptime/alpha)Beta))) 

(2) 

Wherein the Alpha and Beta constants are the customary Weibull constants (i.e., alpha is 
the scale parameter and beta is the shape parameter) known in the art. Equation (2) can be 
20 rearranged to the following form when solved for uptime: 

Uptime = Alpha(-ln (probability of failure)) ( 1 /Beta) (3) 

Equation (2) can be used to calculate uptimes based upon generation of a random 
25 probability, wherein a uniform random real number between zero and one is used for the 
probability of failure. Any log-normal equation known in the art, having parameters 
sigma (i.e., the standard deviation of the log of the downtimes) and mu (i.e., the mean of 
the log of the downtimes), can be used to parameterize a downtime distribution. Purely 
for purposes of discussion herein, the uptime distributions for the failure modes of the 
30 manufacturing system 22, as set forth in Table 1, are parametrized using a Weibull PDF 
and the downtimes are characterized using a log-normal function, wherein the probability 
constant C represents the probabilities among the various downtime distributions. It is 
should be noted that the summation of the constant C values of any single failure mode 

10 
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must equal one, because the downtime from any loss event must fall into one of the 
subject downtime distributions for the failure mode of that loss event. The false start 
probability of each failure mode is represented by a constant P. Illustrative alpha, beta, 
mu, sigma, and C and P constants for each failure mode are set forth below in Table 4. 
As previously discussed, some of the failure modes have more than one downtime 
distribution. For example, a belt driven roller might fail because of failure of a driving 
belt. The roller might be quickly repaired if multiple belts are hung about the driving 
pulleys such that disassembly of the roller structure is not required for the repair and a 
short repair time occurs. However, once the last repair belt has been used, disassembly of 
the roller structure would be required resulting in a longer repair time distribution than if 
a repair belt had been available. Multiple downtime distributions are also useful for 
modeling failure modes having more than one cause, where each cause has a separate 
downtime repair distribution. The number of downtime distributions can be selected 
using one of several techniques known in the art (e.g., segmented regression). 



TABLE 4 



Simulation 


Rate 


Failure 


Failure Mode 


Uptime 


Downtime 


False 


Object 


Constant R 


Mode 




(minutes) 


(minutes) 


Start 




(ft/minute) 


No. 




Alpha/Beta 


C/mu/sigma 


P 


Sheet Roller 


180 


1 


Sheet jam 


36,801/0.8 


0.6/2.6/1.1; 
0.35/5.8/2; 
0.05/43/2 


0.97 






2 


Roll change 


63,400/15 


1/3/1.8 


0.96 






3 


Belt failure 


890,000/14.6 


1/56/2.5 


0.98 


Embosser 


180 


4 


Infeed jam 


345,672/0.6 


1/1.2/0.7 


0.99 


First Cutter 


180 


5 


Blade failure 


200/1.2 


1/180/5.6 


0.97 


Core Roller 


180 


6 


Core misalignment 


947,505/0.4 


1/1.8/0.78 


0.97 


Accumulator 


180 


7 


Accumulator full 


29,023/0.8 


1/3.6/0.98 


1 


Second Cutter 


250 


8 


Blade failure 


3,984/0.93 


1/21/0.9 


0.9 


Polymer roll 


250 


9 


Polymer roll change 


200/1.3 


1/3.9/1.32 


0.9 


Wrapper 


250 


10 


Polymer sheet jam 


2,900/0.87 


1/0.98/0.97 


1 






11 


Polymer sheet 
alignment 


9E+10/0.2 


1/0.798/1.6 


0.98 



Also, listed in Table 4 is a rate constant R which represents a steady state product rate 
value for a simulation object (e.g., 180 feet per minute of sheet feed for the sheet roller). 
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After completing the collection and parameterization of the system loss event data 
at steps 60 and 62 of process 20 of Fig. 1, a computer program modeling or simulating the 
manufacturing system 22 is constructed in accordance with another aspect of the present 
invention. The simulation program is executed or run on any general or special purpose 
5 computer or other digital processing apparatus, such as a desktop computer, a server 
and/or client computer interconnected by a network (for example, either via the Internet 
or an Intranet), a micro-computer, hand-held organizers and other forms of computers and 
computer systems as otherwise known in the art. Exemplary networked computers 1000 
are illustrated schematically in Fig. 5. The computer preferably comprises a logic circuit 

10 (such as central processing unit 1030, microprocessor or other micro controller) capable 
of executing the simulation program. The simulation program, or portions thereof, can be 
provided as a program product, wherein the program product includes a signal bearing 
medium which can be configured to store data and/or machine readable instructions 
which cause the logic circuit to which it is connected to perform the simulation steps 

1 5 discussed hereafter. The signal bearing medium can be provided in the form of an optical 
disk, a magnetic disk, a magnetic hard drive (e.g., reference numeral 1010), a magnetic 
tape, RAM, ROM, or any other magnetic, optical, or other computer readable storage 
medium. Alternatively, the program product can be distributed with the instructions 
contained in other signal bearing media including digital and analog communication links 

20 (e.g., such as a wire or fiber portion of a local area network, a wire or fiber portion of a 
wide area network, a portion of a wireless network, etc.), a carrier wave or propagated 
signal, and other forms of transmission media. The computer also preferably comprises 
one or more input/output peripherals, such as a keyboard, mouse, touch screen, 
microphone, display monitor, printer, etc., which can be interconnected with the logic 

25 circuit via a system bus and adapter (e.g., display adapter 1040). While the step 148 of 
Fig. 1 is preferably computerized, it will be appreciated that the other steps (or portions 
thereof) of the process 20 can be implemented by or in combination with a digital 
processing apparatus. 

Referring to Fig. 6, a preferred simulation program 80 comprises a plurality of 

30 simulation objects 82 and a controller 84, wherein each simulation object 82 represents a 
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physical component of the modeled system (e.g., a plant, a manufacturing line, a machine, 
a component, a queue, etc.). The simulation program can be implemented using any one 
of a number of object oriented or non-object oriented programming languages known in 
the art (e.g., C, C++, EXCEL macros, etc.). In this instance, the simulation objects each 
5 represent a machine 56 of the manufacturing system 22 and comprise instructions and/or 
data which describe the operation of that machine. Thus, there are simulation objects 
illustrated in Fig. 6 for each of the sheet roller 24, the embosser 28, the first cutter 30, the 
core roller 34, the accumulator 38, the second cutter 38, the polymeric roller 42, and the 
wrapper 44. The simulation objects 82 are interconnected with their corresponding 

10 upstream and downstream simulation objects such that a virtual representation of the 
manufacturing system 22 is arranged. The simulation objects 82 preferably send to their 
immediately adjacent upstream and downstream simulation objects the simulated rate 
(e.g., a linear rate in feet/minute of the sheet 26 between simulation object 1 and 
simulation object 2) and the simulated state of the simulation object (e.g., either stopped 

15 or operational). The controller 84 coordinates the simulation time (hereinafter "the 
current simulation time") among the various simulation objects 82. Preferably, the 
coordination is accomplished by each simulation object transmitting its uptime and/or 
downtime to the controller 84 as each simulation object encounters a simulated loss event 
(e.g., a stop). As a simulated loss event is encountered, the controller 84 advances the 

20 current simulation time according to the downtime or uptime calculated at a simulation 
object, as described more fully hereafter. 

Referring to Figs. 6 to 8, an exemplary process 86 implemented in combination 
with a simulation object 82 in accordance with yet another aspect of the present invention 
will now be described. The steps of the process are described both genetically for any 

25 simulation object 82 and by way of example to the sheet roller simulation object, although 
it will be appreciated that these steps can be similarly implemented for the other 
simulation objects 82 of the simulation 80. While the steps of the process 86 are 
described herein with respect to simulation objects, some steps implemented by the 
controller 84 are also incorporated into the process 86 for ease of discussion. It will be 

30 further appreciated that the steps of the process 86 and the arrangement of the simulation 

13 



8070L$&M 



objects and controller can be modified, rearranged, combined, and separated as is known 
in the art without departing from the scope of this invention. 

The process 86 begins at step 88 where the parameters (e.g., constants P/C/R, 
Alpha, Beta, mu, sigma, etc.) for the uptime and downtime distributions, rates, and false 
5 start probabilities of the simulation object 82 are preferably read from an input file or the 
like. Once each simulation object 82 is so initialized, a false start probability array is 
generated at step 89, wherein the array preferably has the structure set forth in Table 5 
below. The first column is the false start probability constant P while the second column 
is the difference of one minus the false start probability P. The third column is the 
10 quotient of one minus the false start probability constant P over the sum of one minus the 
false start probability constant P from column 3. The fourth column is the running 
summation of the third column. 



TABLE 5 



Failure Mode 
Number 


False Start 
Probability P 


1-P 


(l-P)/(sum(l-P)) 


Sum 


1 


0.97 


0.03 


0.3333 


0.3333 


2 


0.96 


0.04 


0.4444 


0.7777 


3 


0.98 


0.02 


0.2222 


1 



15 Next, a target false start probability is summed for the competing cause failure 

modes of the simulation object according to the following equation at step 90: 

Target false start probability = 1- S (1- false start probabilities P) (4) 

20 After determining the target false start probability at step 90, a uniform random real 
number between zero and one is generated at step 94 and this uniform random real 
number is used to calculate an uptime for each cumulative cause failure mode for the 
subject simulation object 82 based upon equation (2) above. For example, the sheet roller 
24 has the roller change and belt failure cumulative cause failure modes, as listed in 

25 Tables 1 and 4. The illustrative alpha and beta values for failure mode number 2 listed in 
Table 4 are 63,400 and 15, respectively. For a uniform random real number value of 0.54, 
this cumulative cause failure mode would have an uptime of 63,346 minutes, respectively, 
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at step 94. These uptime values are preferably stored in an array or other data structure 
which is updated as the simulation 80 is executed. 

Steps 96 to 106 are next executed to calculate a competing cause failure mode 
false start error. At step 96, a loss event counter storing the total number of simulated 

5 loss events (i.e., both competing cause and cumulative cause failure mode loss events) is 
incremented as each loss event is encountered. If the previous loss event for the 
simulation object was due to a competing cause failure mode and was a false start as 
indicated at step 98, step 100 is executed, wherein a counter summing the number of 
competing cause failure mode false starts is incremented by one. Otherwise directly 

10 execute step 102, wherein the simulation's false start event probability is calculated 
according to the following equation: 

Simulation false start prob. = T otal number of CCFM false start events (5) 

Total number of simulation loss events 

15 

The total number of simulation loss events in the denominator of equation (5) was 
previously calculated at step 96, and the total number of competing cause failure mode 
false start events in the numerator of equation (5) was previously calculated at step 100. 
20 The error between the simulation's false start probability and the target false start 
probability is next calculated according to the following equation at step 106: 

Error = Simulation false start probability - Target false start probability 

25 Once the error has been calculated, a corrected simulation false start probability is 
calculated at step 1 06 based upon the calculated error. Any error correction methodology 
known in the art can be employed to calculate the corrected simulation false start 
probability value (e.g., a binary search). While an error correction is described herein for 
simplicity with respect to false starts following a competing cause failure mode, it will be 

30 appreciated that an error correction could be similarly implemented for false starts 
following a cumulative cause failure mode. 
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At step 110, it is determined whether the previous loss event was due to a 
competing cause failure mode, or whether it was due to a cumulative cause failure mode. 
If it is due to competing cause failure mode, steps 112 and 114 are used to determine if 
there is a false start following the downtime (i.e., following the repair time) for the 

5 previous loss event while accounting for the corrected simulation false start probability 
determined at step 106. If the previous loss event is due to a cumulative cause failure 
mode, then steps 111 and 113 (Fig. 8) are executed. At step 1 1 1, a uniform random real 
number between zero and one is generated. If the uniform random real number is less 
than or equal to the false start probability P of the previous cumulative cause failure 

10 mode, as described at step 113, then there is a false start following this loss event. The 
downtime for this false start event is next generated using process 116 (Fig. 9), as 
discussed more fully hereafter; otherwise there is no cumulative cause false start 
following the previous cumulative cause loss event and step 112 is next executed to 
determine if there is a competing cause false start event beginning with step 112, as also 

15 discussed more fully hereafter. 

At step 110, if the previous stop was for a competing cause loss event (or if step 
110 is executed for the first time at startup of the simulation), steps 112 and 114 are 
performed. A uniform random real number between zero and one is generated at step 112 
and if this uniform random real number is greater than or equal to the corrected 

20 simulation false start probability determined at step 106, then there is no false start 
following the previous loss event, then parallel paths A and B (Fig. 8) are performed to 
determine the next minimum uptime for the subject simulation object 82. If the uniform 
random real number generated at step 114 is greater than the corrected false start 
probability determined at step 106 (i.e., there is a false start event following the previous 

25 loss event), then a new uniform random real number is generated at step 115 and this 
number is compared to the fourth column of the array generated at step 89 to determine 
which failure mode of the subject machine step has caused the false start. For example 
with reference to Table 5, if the uniform random real number generated at step 1 1 5 is 
0.85, then failure mode number 3 would be used to determine the downtime for the false 

30 start event because 0.85 is greater than 0.7777 but less than 1. After determining which 
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failure mode causes the false start event, the downtime process 116 (Fig. 9) is next 
executed to determine the downtime for the false start event. 

Referring to Fig. 8, parallel path A calculates the uptimes for each competing 
cause failure mode of the subject simulation object while parallel path B generates the 
5 uptimes for each cumulative cause failure mode of the subject simulation object. 
Beginning at step 1 1 8 for path A, a uniform random real number between zero and one is 
separately generated for each competing cause failure mode of the simulation object. 
Each uniform random real number is applied to equation (3), wherein equation (3) 
incorporates the alpha and beta for the failure mode corresponding with the generated 

10 uniform random real number. For example, Table 6 lists the uniform random real number 
generated for each competing cause failure mode of the sheet roller 24. Using the alpha 
and beta for each failure mode, an uptime is calculated based upon the generated uniform 
random real number for a competing cause failure mode. The smallest competing cause 
failure mode simulation time is next selected at step 120 from among all the competing 

15 cause failure mode simulation times for the subject simulation object and this value is 
forwarded to step 132. 



TABLE 6 



Failure 
Mode 


Alpha/beta 


Random 
No. 


Current 
Simulation 
Time 


Uptime 
(minutes) 


Competing cause failure 
mode simulation time 


Sheet jam 


36,801/0.8 


0.14 


10 


84,308 


84,318 



Referring now to parallel path B of Fig. 8, the smallest uptime among the 
20 cumulative cause failure modes for the subject simulation object 82 is determined at steps 
122 to 130. At step 122, for each cumulative cause failure mode of the subject simulation 
object 82, the system uptime is subtracted from the cumulative cause failure mode uptime 
of each cumulative cause failure mode and these new values are stored in the array. At 
step 124, it is determined whether the cause of the previous loss event is a cumulative 
25 cause. If yes (e.g., a cumulative cause failure mode has caused a stop), then a new 
cumulative cause failure mode uptime is calculated for that cumulative cause failure mode 
by generating a uniform random real number between zero and one at step 126; otherwise, 
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step 130 is executed as discussed below. Using equation (3) and the alpha and beta for 
the subject cumulative cause failure mode, the new uptime is calculated which replaces 
the expired cumulative cause failure mode uptime in the array. The smallest cumulative 
cause failure mode uptime is next selected from among all the cumulative cause failure 
5 mode uptimes in the array for the subject simulation object 82 at step 130. Next, the 
smallest uptime value between step 130 and step 120 is selected at step 132 to determine 
which competing cause failure mode or cumulative cause failure mode will cause the next 
loss event for the subject simulation object 82. 

After determining which failure mode causes the next loss event for the simulation 

10 object, this smallest failure mode uptime is sent to the controller 84 at step 133. At step 
135, the controller 84 selects the smallest failure mode uptime from among all the 
simulation objects (the "system uptime"), including from any accumulator product level 
loss events (discussed more fully hereafter as another preferred embodiment), and adds 
that uptime to the current simulation time, thereby advancing the simulation to the 

1 5 beginning of the loss event. The system uptime is returned to each simulation step. At 
step 139, each simulation object determines whether its failure mode uptime determined 
at step 132 matches the returned system uptime. If it does, then execution passes to the 
downtime process 116 which is performed beginning with step 134 (i.e., that simulation 
object is the one encountering the loss event), as shown in Fig. 9; otherwise, steps 141 

20 and 143 are performed wherein the simulation object enters a "wait or idle" state at step 
141, wherein each simulation object awaits input from its adjacent simulation objects 
and/or the controller. Referring to Fig. 9 and the downtime process 116, updated rate and 
state information is sent to each adjacent simulation object at step 134, which generally 
for a loss event is a rate of zero and a state indicator that the subject simulation object 82 

25 has stopped. The state indicators are merely used to coordinate events among the 
simulation objects. A uniform random real number between zero and one is generated at 
step 136, and this number is used at step 138 to select which downtime distribution will 
be used to calculate the downtime for the particular loss event by comparing each 
downtime constant C to the uniform random real number. The constants C are 

30 sequentially summed and the uniform random real number is compared to the range in 
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between the summations to determine which downtime distribution to use. For example 
with reference to Tables 4, failure mode 1, and Table 7, if the uniform random real 
number generated at step 136 is 0.85, then the downtime distribution # 2 from Table 4 
would be used to determine the downtime because 0.85 is greater than 0.6 but less than 
5 0.95. 



TABLE 7 



Downtime 
Distribution No. 


Constant C 


Summation of 
constant C values 


Range between 
summations 


Random No. 


1 


0.6 


0.6 


0< random number < 0.6 




2 


0.35 


0.95 


0.6<random number<0.95 


0.85 


3 


0.05 


1.0 


0.95< random number 





At step 140, a new uniform random real number between zero and one is 
generated. This uniform random real number is used with equation (4) to determine the 

10 downtime for the loss event at step 142 using the downtime distribution parameters for 
the downtime distribution determined at step 136. The downtime calculated at step 142 
(the system downtime) is sent to the controller 84 at step 144. The controller adds the 
system downtime to the current simulation time (i.e., the simulation time is now advanced 
to the end to the repair period) and forwards the current simulation time to each of the 

15 simulation objects at step 144 (this information being received by the simulation objects 
at step 143). The rate and state for the subject simulation object 82 which caused the loss 
event are next updated and sent to the adjacent simulation objects at step 146, wherein 
generally the rate is set to the rate constant (e.g., 180 ft/minute for the sheet roller 24) and 
the state is updated to running. Execution, for all simulation objects 82 then resumes at 

20 point E of Fig. 6, wherein the false start probability error correction begins at step 96. 

An alternate preferred embodiment of the present invention incorporates the 
process 152 illustrated in Fig. 10 for simulation objects modeling storage or buffer 
components, such as the accumulator 38. Each accumulator has an empty product level 
value (i.e., the value at which there is no product in the accumulator), a current product 

25 level value (i.e., the current value for the amount of product in the accumulator), and a 
full product level value (i.e., the value at which the accumulator is completely filled with 
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product and can accept no more). Beginning at step 154, the difference between the rate 
of product into the accumulator and the rate of product out of the accumulator is 
calculated. By way of example with respect to the simulation 80, the rate of product into 
the accumulator 38 of the manufacturing 22 is the difference between the rate constant for 

5 product delivered from the simulation object 82 representing the core roller 82 and the 
rate constant for the product received by the simulation object 82 representing the second 
cutter 40. If the difference is greater than zero at step 156, then step 158 is performed, 
wherein it is determined how long it is before the accumulator 38 is full by subtracting the 
full product level value from the current product level value and dividing this difference 

1 0 by the rate difference calculated at step 1 54. This value is then summed with the current 
simulation time received by the simulation object 82 for the accumulator at step 143 of 
Fig. 6 and sent to step 132. This accumulator level simulation time value (which 
represents a loss event) is then compared to the smallest competing cause and cumulative 
cause simulation times of step 132. Referring to step 160, if the rate difference is less 

15 than zero then step 162 is executed, whereat it is determined how long it is before the 
accumulator 38 is empty by subtracting the current product level value from the empty 
product level value and dividing this difference by the rate difference calculated at step 
154. This value is then summed with the current simulation time received by the 
simulation object 82 for the accumulator at step 143 of Fig. 6 and sent to step 132. This 

20 accumulator level simulation time value (which represents a loss event) is then compared 
to the smallest competing cause and cumulative cause simulation times of step 132. 
Otherwise step 164 is performed (i.e., the rate difference equals zero and the rate in equals 
the rate out), wherein the product level within the accumulator is constant and an 
accumulator level simulation time value of infinity is sent to step 132. 

25 Referring to Table 8, six exemplary time steps (e.g., T=0 minutes, T =143 

minutes, T=322 minutes, T=327 minutes, T=330 minutes, and T=350 minutes) for a 
simplified simulation in accordance with the present invention of the manufacturing 
system 22 will now be described. Each time step incorporates one or more steps of the 
previously described simulation processes and is intended to generally demonstrate the 

30 concepts described herein for further clarity. At the first time step (i.e., at T=0 minutes 
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or the beginning of the simulation after initialization of the simulation), a uniform random 
real number is generated at the fourth column to determine whether there is a false start 
for the competing cause failure modes under the simplified assumption that there is not a 
cumulative cause false start at startup. The same single random real number is used for 

5 each competing cause failure mode of a machine object (i.e., steps 1 12 and 114 of Fig. 7) 
step and a separate random real number is generated for each cumulative cause failure 
mode when appropriate (i.e., when step 1 1 1 of Fig. 7 is executed). If there is a false start 
as determined at steps 113 or 114 of Fig. 7, a "yes" is indicated in fifth column. Since 
there is no false start for the first time step, a uniform random real number is generated for 

10 each failure mode and this uniform random real number is then used to generate an 
uptime. The smallest uptime for each simulation object is selected at the eighth column 
and then the smallest uptime for the system is selected at the ninth column. The system 
uptime is added to the current simulation time at the second time step (i.e., time T= 0+143 
minutes). A uniform random real number is generated at the tenth column; however since 

15 there is only one downtime distribution, it is not does not affect which downtime 
distribution is selected. At the eleventh column, a random number is generated, which is 
then used to determine the downtime for failure mode number five. At the third time step, 
the downtime from the twelfth column is added to the simulation time (i.e., time T = 143 
+ 179). Uniform random real numbers are again generated at the fourth column for the 

20 false start probabilities (except the cumulative cause failure modes since the prior loss 
event was not due to a cumulative cause failure mode), and there are no false starts. New 
uptimes are generated for the competing cause failure modes only at seventh column 
based upon the uniform random real numbers listed in the sixth column. The uptimes for 
the cumulative cause failure modes at the seventh column are equal to the uptime from 

25 the first time step (i.e., at time T = 0) minus the system uptime calculated at the first time 
step. A new system uptime is selected at the tenth column (e.g., 5 minutes). Following 
the downtime calculation of fourth time step (i.e., at time T = 327 minutes), a false start 
occurs at the fifth time step (i.e., at time T = 330 minutes), whereat a downtime is 
calculated for the failure mode eight which caused the false start event. This is calculated 

30 according to steps 134 to 146 of Fig. 9. At the sixth time step (i.e., time T = 350 
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minutes), the downtime for the false start event calculated at the fifth time step is added to 
the current simulation time at the first column of the fifth time. The uptime for each 
competing cause failure mode is generated at the seventh column based upon the uniform 
random real numbers of the sixth column. Since there was no system uptime between the 
fourth and sixth time steps, the uptimes for the cumulative cause failure modes remain the 
same as the uptimes of the third time step. 



TABLE 8 



Simulation 
time 


Simulation object 


Failure 
mode 
no. 


Random no. 
for false start 
probability 


False 
start? 


Random 
no. for 
failure 
mode 
"uptime 


Pailure mode 
uptime 


smallest failure 
mode uptime of 
object 


nallest 
ailure 
mode 
ptime 

of 
ystem 


Random 
no. for 
owntime 
stribution 
selection 


tandom 
no. for 
>wntime 


owntime 


T=0 


Sheet Roll 


1 


0.23 


No 


0.14 


84306 
















2 


none 


N/A 


0.54 


63346 


63346 














3 


none 


N/A 


0.83 


598295 














Embosser 


4 


0.42 


No 


0.63 


90404 


90404 












First Cutter 


5 


0.08 


No 


0.51 


143 


143 


143 










Core Roller 


6 


0.62 


No 


0.18 


3575669 


3575669 












Accumulator 


7 


0.057 


No 


0.47 


20011 


20011 












Second cutter 


8 


0.61 


No 


0.37 


3906 


3906 












Polymer 


9 


none 


N/A 


0.5 


148 


148 












Wrapper 


10 


0.8 


No 


0.5 


1888 


1888 














11 


0.8 




0.79 


54117352 


54117352 










T=143 


Sheet Roll 


1 
























2 
























3 






















Embosser 


4 






















First Cutter 


5 














0.38 


0.41 


179 




Core Roller 


6 






















Accumulator 


7 






















Second cutter 


8 






















Polymer 


9 






















Wrapper 


10 
























11 




















T=322 


Sheet Roll 


1 


0.38 


No 


0.71 


9550 
















2 


none 


N/A 


None 


63203 


63346 














3 


none 


N/A 


None 


698152 














Embosser 


4 


0.73 


No 


0.48 


198272 


90404 












First Cutter 


5 


0.34 


No 


0.1 


400 


400 












Core Roller 


6 


0.82 


No 


0.32 


1267862 


1267862 












Accumulator 


7 


0.16 


No 


0.57 


14035 


14035 
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Second cutter 


8 


0.28 


No 


0.14 


8207 


8207 












Polymer 


9 


none 


N/A 


None 


5 


5 












Wrapper 


10 


0.65 


No 


0.4 


2622 


2622 














11 


0.65 


No 


0.16 


7.49E+09 


7.49E+10 










T=327 


Sheet Roll 


1 
























2 
























3 






















Embosser 


4 






















First Cutter 


5 






















Core Roller 


6 






















Accumulator 


7 






















Second cutter 


8 






















Polymer 


9 














0.75 


0.29 


3 




Wrapper 


10 
























11 




















T=330 


Sheet Roll 


1 


0.61 


No 




















2 


none 


N/A 




















3 


none 


N/A 


















Embosser 


4 


0.38 


No 


















First Cutter 


5 


0.32 


No 


















Core Roller 


6 


0.78 


No 


















Accumulator 


7 


0.487 


No 


















Second cutter 


8 


0.99 


Yes 










0.45 


0.41 
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Polymer 


9 


0.31 


No 


















Wrapper 


10 


0.7 


No 




















11 


0.7 


No 
















T=350 


Sheet Roll 


1 


0.7 


No 


0.02 


63197 


3197 














2 


none 


N/A 


None 


63203 
















3 


none 


N/A 


None 


698152 














Embosser 


4 


0.83 


No 


0.86 


10021 


10021 












First Cutter 


5 


0.81 


No 


0.21 


288 


288 


288 


0.53 


0.65 


182 




Core Roller 


6 


0.05 


No 


0.61 


152322 


152322 












Accumulator 


7 


0.45 


No 


0.92 


1300 


1300 












Second cutter 


8 


0.67 


No 


0.98 


378 


378 












Polymer 


9 


none 


N/A 


None 


488 


488 












Wrapper 


10 


0.56 


No 


0.77 


625 


625 














11 


0.56 


No 


0.46 


2.54E+10 


2.54E+10 











Referring again to Fig. 1, after the simulation 80 is built at step 148 as described 
above, the simulation 80 can be used to analyze the manufacturing system 22 and 
implement changes in the manufacturing system 22 to improve its reliability, as described 
at step 150. For instance, the size of the accumulator 38 might be changed, the rate 
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constants for a machine might be adjusted, the arrangement of machines in the 
manufacturing system might be changed, the value of eliminating a failure mode, etc. In 
addition, the system data can be used to model a second, yet to be built, manufacturing 
system having similar system objects, perhaps arranged in a different manner, in order to 

5 predict and optimize the reliability of the second system prior to building it. 

Set forth below are examples of availability, MTBF, and MTTR for an existing 
manufacturing system and the same values from a simulation made in accordance with 
the present invention, wherein the simulation modeled cumulative and competing failure 
modes with beta values greater and less than one and false start probabilities. Preferably, 

10 the percent error (i.e., the ratio of the actual value minus the simulated value over the 
actual value X 100) is less than about 3%, and more preferably less about 2%. Most 
preferably, the error is less than about 1%. As will be appreciated, other values, such as 
the number of loss events, the downtime by loss event, can be outputted from a simulation 
made in accordance with the present invention. The amount of the error can be further 

15 decreased by, for example, increasing the number of downtime distributions which are 
simulated. 



Example 1 





Actual system 


Simulation of actual system 


Percent Error 


Availability 


0.8069 


0.8067 


-0.1 


MTBF 


14.30 


14.29 


-0.1 


MTTR 


3.39 


3.42 


1 
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Example 2 





Actual system 


Simulation of actual 
system 


Percent Error 


Availability 


0.8616 


0.8625 


0.1 


MTBF 


17.40 


17.86 


2.6 


MTTR 


2.80 


2.85 


1.7 



Example 3 



Failure 


Actual 


Simulated No. 


Actual 


Simulated 


Percent Error 


Percent 


Mode 


No. 


of Events 


Downtime 


Downtime 


In No. of Events 


Error in 




Events 










Downtime 
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1 


4192 


4271 


7355 


7116 


-1.89 


3.24 


2 


3012 


3056 


18018 


18514 


-1.48 


-2.75 


3 


752 


782 


1178 


1130 


-3.98 


3.99 


4 


387 


400 


48 


489 


-3.44 


-1.06 


5 


143 


145 


96 


100 


-1.63 


-3.92 


6 


62 


61 


82 


84 


0.53 


-2.54 


7 


39 


39 


136 


142 


-1.53 


-3.90 


8 


99 


100 


108 


105 


-1.21 


2.76 



The foregoing description of the preferred embodiments of the invention have 
been presented for purposes of illustration and description. It is not intended to be 
exhaustive or to limit the invention to the precise form disclosed. Modifications or 
variations are possible and contemplated in light of the above teachings by those skilled 
in the art, and the embodiments discussed were chosen and described in order to best 
illustrate the principles of the invention and its practical application. It is intended that 
the scope of the invention be defined by the claims appended hereto. 
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